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Abstract 

A mechanism generating a A-type behavior in the specific heat of a Bose gas 
near the critical temperature T c is discussed. It is shown to work for a general 
class of quasiparticle spectra, and especially, for the Bogoliubov's model of a weakly 
interacting Bose gas, where a temperature-dependent quasi-particle spectrum is 
obtained. 



1 Introduction 



The first experimental measurement of the low temperature specific heat of He 4 was 
made by Keesom in 1927-8 [|TJ who became impressed by its A shape and by its singular 
behavior at the critical temperature where superfluid behavior started to manifest. In 
his book on Statistical Mechanics, Feynman refers to the A point behavior as one of 
the unsolved problems of superfluidity theory, and expresses his view that perhaps part 
of the explanation of the lambda transition involves Bose condensation since the c v (T) 
-curve for an ideal Bose gas shows a peaked form at the critical temperature, although it 
has not a divergent behavior. We are not aware about any satisfactory model explaining 
this behavior, except the one presented in the paper by D.M.Ceperley || in which it is 
simulated the boson system by means of the Monte Carlo techniques using path integrals. 
In that paper, a good agreement is obtained among the simulation and the experimental 
measurements of the He 4 specific heat. 

In our present letter we want to discuss a mechanism for generating a A-type behavior 
in the specific heat of a Bose gas. We consider at first some models for the phase transition 
of symmetry restoration in temperature relativistic systems with spontaneous symmetry 
breaking, which exhibit a A-type behavior for c v , and then we study the similar problem 
in the non-relativistic case. We show that for the weakly interacting Bose gas one can 
find such a divergent behavior of the specific heat as being produced by the condensate. 
Obviously, such a system is far from being a satisfactory model for the superfluid, in 
which the interactions cannot be considered weak at all. However, as the mechanism is 
essentially an infrared (long wavelength) effect, it may provide insights in understanding 
what happen in a more exact model of superfluidity. Also, the behavior of the specific 
heat for the weakly interacting Bose gas may become interesting in connection with the 
recent experimentally observed Bose condensate, although such problem would require to 
consider the particles in an external field, provided by the magnetic trap. 



neous symmetry breaking (SSB). For the latter, one can find that a singular behavior of 
the specific heat occurs at the symmetry restoration temperature, due to the temperature- 
dependent mass, proportional to the symmetry-breaking parameter. We prove that a sim- 
ilar phenomenon occurs for boson systems whose quasiparticle spectra have some specific 
non- vanishing temperature-dependent spectra in the infrared limit. 

We will take either the simple model of the scalar field or the Abelian Higgs model 
||, [Q. One can write the effective potential as 

v(0 = ^-^f + v(T,o, (i) 



where V(£) is the sum of tadpole diagrams. The spectrum is e(p)i = c^j {p 2 + M 2 c 2 , 
where i = 1 for the scalar and i = 1,2,3,4 for the multiplet in the scalar-vector Higgs 
model where, Mi = A£, and M2^,a = g£,/2 and g is the scalar-vector coupling constant. 
In the high temperature limit is V(T, £) = <xT 2 £ 2 /2, the extremum of V(£) leading for 



T < T c to a dependence M(T)j = Ki\jT 2 — T 2 , T c being the critical temperature for 
symmetry restoration. The infrared contribution to the thermodynamic potential in the 
limit T ^> M(T)i is Afl ~ i2~d — • O ne obtains then a specific heat 

Ac, = -T%^ ~ Vc\ 3 f (2) 

dT 2 ^ l 2ixWjff~T 2 KJ 

which diverges for T — > T c , showing a A-type behavior. 

In general, we may argue that any quasi-particle spectrum of the Bose gas whose 



infrared limit behaves as £:(p) = ay/5' 2 + p 2 , where 5 = 6y 1 — (T/T c ) 7 , and 7 > 1, leads 
to a A-type behavior of c v . One can obtain that result from the infrared term in the 
calculations made in [Q. But also, by taking T close enough to T c to make 5 arbitrary 
small, one can choose some small momentum rj 3> 5. Then the integral over p in the 
thermodynamic potential Q = J ln[l — exp(— (3e(p))} can be decomposed in the 
sum Jq + J°°, from which the infrared contribution as AQ = 5 3 T/12nh 3 is obtained, 
leading to the A-type behavior for c v = —Td 2 Q/dT 2 . 



Obviously, T < T c , where T c is the critical temperature for condensation and to describe 
the system we use a version of the Bogoliubov's model |J valid near T c , which leads to a 
spectrum having the abovementioned infrared limit. 

2 The Bogoliubov's model near the critical point. 

We shall assume an weakly interacting Bose gas under the approximate conditions a/A 
1 and na 3 <C 1 where a is the scattering length of the two body interaction, A = h/\ / 2iTmT 
the thermal wavelength of the particles and n the particle density in the system considered. 
We will start from the quantized field Hamiltonian expanded in terms of the quantities 
dp, a p (p 7^ 0) as in the usual Bogoliubov's procedure 0. In what follows we refer briefly 
to the details of this model by following Pathria f7j. 

Our case differs from Bogoliubov's one by two facts :1-) as our temperature interval 
is close to the critical temperature, we cannot consider in all cases that the number of 
particles in states with p ^ is always much more smaller than n , number of particles 
in the condensate. 2-)The usual approximation n ~ N cannot be made. 

Thus we start also from the Hamiltonian of the quantized field for spinless bosons, we 
assume momentum conservation in the interactions and name U(r) the repulsive potential 
of the two body interaction, Upl'p% = y J J d 3 re^~U(r) 

The scattering length can be written, then, as a = 

z5&£#,m, and as we are con- 
sidering a repulsive potential always a > 0. For p = (zero momentum transfer in the 
colision) U D = J d 3 rU(r). 

In the temperature interval we are considering the momenta are very small and we 
can assume that the momentum transfer in each collision is almost zero p ~ 0. For that 
reason it is possible to express approximately the matrix elements by using U a . Then the 
Hamiltonian is, 

p 2 U 

^ = E^r a p a p + ^7E<' 1 a p' 2 a P2«Pi- ( 3 ) 



a a ~ n 2 Q which allows us to take each one of these creation and annihilation operators 
as c-numbers proportional to ■JrT - Then the second term in the Hamiltonian can be 
expanded, its lowest-order term being a^a^a a = n 2 — n and the lowest-order ground 
state energy is E Q = ^[n 2 — n ). The next order terms to be included leads to the term 
{ a t a t + a p a -p + 4ai"a_ p }n and our final Hamiltonian is 

J2i a t a - P + a P a -p + 4a p a -p) • 

The usual Bogoliubov's procedure is, by considering that n ~ N, (which ceases to 
be valid near the transition temperature) to replace N by no in the in the third term in 
(HI) (this means to neglect one term of fourth order in the quantities ap~,a p ), and after 
introducing the relation a P a p + n o — N in the first, leads to the cancellation of a term 

p^O 

^a p a± p in the last bracket of (g). 

Thus, our assumptions lead simply to change the coefficient of the last term in curly 
brackets from 2 to 4. The next step is to make the usual Bogoliubov transformation 



a p — (bp — ctp&i p )/ y/l — a p 2 and a+ = (6+ — a p b_ p )/Jl — a p and as a result of it we get 
the diagonalized Hamiltonian, 



H = E + J2<p)b;b p , (5) 

p^O 

where E = ^[n 2 Q - n a ] - U ^r E " P and a p = ^[4%^ + ^ - e(p)], and if we name 

p^O 

K = ^ then 



is the spectrum of the new Bose quasiparticles representing the elementary excitations of 



In Bogoliubov's spectrum the term 12K 2 is not present. In our model the long wave- 
length limit is obviously not linear in p. As K <C kT is very small (for no — 10 16 , K ~ 
1CT 12 eV) it can be usually neglected, but it is able to produce the typical A-type di- 
vergent behavior of the specific heat near T c . As lim p ^o^(p) — 2y/3K, the parameter 
K formally behaves as the analog of a rest energy in relativistic dynamics. This "rest 
energy" has the remarkable property that it decreases with temperature and goes to zero 
for T — > T c . Actually, K is proportional to the symmetry breaking parameter no, which 
is the condensate, and condensation means some (first order) gauge symmetry breaking 
Q. Thus K arises as some sort of non-relativistic analog of the temperature scalar and 
abelian Higgs models mentioned above. An explicit computation of c v for this case shows 
a divergent behavior close to the critical temperature, as it is easily checked by doing the 
temperature expansion of the thermodynamic potential, which we will do by taking the 
exact spectrum (H). 

The divergent behavior of the specific heat is related to the fact that at the critical 
temperature a macroscopic number of particles suddenly falls to the ground state. Due 
to the interaction term, the ground state energy is different from zero, 



2irh 2 anl I 16 /6a 3 n„ 
E n = — — < 1 



mV 15 V 



(7) 



where E, F are the usual elliptic integrals. Due to interactions there is also a macroscopic 
number of particles in states very close to the condensate. 

3 Quasi-particle thermodynamic potential 

The thermodynamic potential of the quasi-particles whose spectrum was obtained in the 
last section is f2 = J ln[l — exp(— (3e(p))} and (3 = We have taken the chemical 
potential ji = 0. We are going to do an asymptotic expansion of this potential close to T c 
for T < T c , taking into account that for T — > T c , n — > 0. 



n 



2(2mf/ 2 V 



J dxx 2 ln[l — e" 



v / i2l 5 +8M?T?l 



(8) 



where M = Kf3 <C 1 and we are going to do our expansion in terms of it, in the same 
way done in the case of the effective potential in the abelian Higgs model ||. 
Thus, 



Q = fi(M = 0) + 



on 



dM 



M + R(M), 



(9) 



M=0 



where we stop our expansion in the first-order term and R(M) is certain function of M 
that we are going to find out as 



R{M) = fi(Af) - Q(M = 0) 



dM 



M, 



M=0 



if the first derivative of this expression is calculated we obtain 



(10) 



dR dQ 



M=0 



dM dM dM 
and this is the expression to be used for computing R(M). 

After some calculations that are shown in the appendix the thermodynamic potential 
is obtained as 



n{M) = ¥§^{t t{ 1 K ^ (12) 



3 



3 



const. 

12V7T 3 



The most interesting term for us is the third one, containing the infrared contribution. 
Such behavior can be obtained also by cutting the interval of integration in (ffi) in the 
form / °° = Jq + C° by taking some small rj 3> M. Then in the first integral we can 
expand the exponential in the denominator, leading to a linear in T term, which contains 
the third term of (|13[). 



4 Internal energy and specific heat 

As K — > near T c , in what follows we are taking into account only the terms up to K 3 / 2 
Then for U = - T§p we obtain 



and for the specific heat c v = |^ the expression 



(2m) 3 / 2 /5 

c v 



U {^ r O c 4 )A;5/2T3/2 - 3r ^K(^ 3/2Tl/2 ^ 



h%2 



+4fcT[2arf/ 2 _ 3r( 3 )c( 3 )A;1 / 2Tl/2] ^ 
-4fcT 2 [r(^)C(^ 1 / 2 T^ - avK^f^ 

+a2nkT 2 {^j V 1 / 2 

and it is the last term in (gj) the one giving the divergent behavior of c„ for T — > T c _. 
The first term gives exactly the contribution of the ideal boson gas. 
The contribution of the last term can be expressed in terms of n Q as 



M™\3/2 ( 



By taking n = N 1 — (|r) 7 , (7 > 0, for the ideal gas 7 = |), one obtains that for 
T^T C _, Ac v ~ ... N ^XL 



T 2 7-2 C y/l-(T/T c p ' 

The behavior of c v for T > T c would require a separate investigation. For exactly 



T = T c , as argued in a previous paper |TIJ, for the ideal gas not in the thermodynamic 
limit, the number of particles in the condensate is not zero, and for T ~ T + we may expect 
that the densities in states close to the ground state are large, but rapidly decreasing with 
increasing T. In that region, the chemical potential cannot be taken as zero. Thus, we 
expect that a mechanism similar to the one described for T < T c may perhaps leads to 
a divergent behavior in that region for T — > T c . We shall take, however, the free gas 
behavior in our present approximation. We will write also 

n = Nf(T). (16) 

The specific form of f(T) can be taken as the ideal gas one since usually the quantity 
«p « 1 and thus the average density of quasiparticles n p is of the same order of the 
particles not in the ground state N p (since N p = n p (l + a p /l — a p ) + a p /l — a p ), p"T|j , 
and these are proportional to T 3//2 . We take, thus, 

/ T \ 3//2 ~ 

j(T) =[i-y j, in) 

along with helium parameters, the divergent curve c v (T) for T = T c shown in figure 1 is 
obtained. In this case the scattering length has been chosen as the radius of a hard-sphere 
as a = 2.14A (Kalos,Levesque, and Verlet)[|T^]. We want to remark that such a divergent 
behavior cannot be obtained by using the usual Bogoliubov's spectrum |13| . 

5 Conclusions 



From our results we conclude that some quasi-particle boson spectra having nonvanish- 

incr tpmnprntnrp-rlpnpnrlpnt infrnrprl limit for snmp snppifip HpnpnHpnpps of this limit nn 



the weakly interacting gas it is seen that if one does not take the approximation n Q ~ N, 
but use instead no(T), one obtains a quasiparticle spectrum which is nonvanishing in the 
long wavelength limit, lim p ^ ^(p) = Vl2K where K is proportional to the condensate 
density n . The nonvanishing value of K is a manifestation of the symmetry breaking, 
and for n Q = the spectrum is reduced to the free particle one, where lim P ^<> e(p) = 0. 

1/2 

The infrared contribution to the internal energy contains a term proportional to n . For 
a dependence of n on temperature of form n = N[l — (T/T c ) 7 ], 7 > 0, a A-type divergent 
behavior of the specific heat c v is obtained as T — > T c . 
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7 Appendix 

7.1 Lowest and first order terms of Q 

If f2 is evaluated for M = the expression for the thermodynamic potential is 



tt(M = 0) 



-2(2mf/ 2 V 5 5 
3p 5 / 2 h 3 (2n) 2 l 2 M 2 



(18) 



Then 




dM (3 5 / 2 h 3 {27i) 2 J V12M 2 + 8Mx 2 + x 4 [exp(Vl2M 2 + 8Mx 2 + x 4 ) - 1] ' 



on 



DM 



M: 



4(2m) 3 /V 3 3 
=0 ^l 2 h\2n) 2 l 2 M 2 J ' 



(19) 



7.2 The R(M) function 

If we substitute in (O) the expressions (113) , (|T^) and the Matsubara sum 



we obtain 



£:[exp(£:) — 1] 



E 



e 2 + Ann 2 2e ' 



(20) 



dR 
dM 



8(2m) 3 /y 
P 5 / 2 h 3 {2Ti) 2 

oo 

+2E 



n=l 

oo 



dxx 2 [3M + a; 2 ] 
12M 2 + 8Mx 2 + x A 

dxx 2 [3M + x 2 } 



dxx 4 



:r' J 



Aim 2 + 12M 2 + 8M:r 2 + x 4 



47rn 2 + x 4 



rfxx 2 [3M + x 2 ] 
V12M 2 + 8Mx 2 + x 4 



dxx 4 



x- 



(21) 



Each one of the terms between brackets are calculated separately, thus the first one is 



-an 



M 1/2 , (22) 



where a = (f + ^ 1 )(v / 6 - y/2). 

In the second term as M is very small in our range of temperatures we can take the 
limiting case of n large and expand it in a power series of — so 



SMJtt 1 M 2 
...i = — ^ + -—-^ + .... 23 



The last one, after integrating is, 

[...], = —ML + f-M^ { 7 £<* yf ) - 2F(* /|)} , (24) 

where the divergent integrals have been been regularized by introducing a cut-off L large 
enough. 

The first term of when the sum over n is done is a divergent series. Thus we take 
an upper cut-off term at n — N(L) such as to cancel both divergent terms. At the end 
we get 



OR 8(2m)^V / tt 1/2 " -3M^ 1 Q(M 2 ) 3 

2 



-\\-ml + f-M^ m l /|) - ml /f: 



))] , 



then the sum can be approximated by 



N 1 ^ rfn 



in this way by taking L = -^LiV 1 / 2 , the otherwise divergent terms cancel each other and 
after the integration in M has been performed the value of the R(M) function is 



R(M) = * (-^M3/ 2 

1 ; 8 5 / 2 h 3 2tt 2 I 3 



-\^ Mii ^i§-<4 )] (27> 

A/ 3 0(M 3 ) 1 



where const does not depend on M. 
R pfprpn r(±<z 



[2] R.P. Feynman, Statistical Mechanics Benjamin, Reading, Mass. (1982). 
[3] D.M. Ceperley, Rev. Mod. Phys. 67 (1995) 2. 
[4] L. Dolan and R. Jackiw, Phys. Rev. D 9 (1974) 3320. 
[5] A.D. Linde, Rep. Prog. Phys 42 (1979) 
[6] N. N. Bogoliubov, J. Phys. USSR 11, (1947) 23. 
[7] R. K. Pathria, Statistical Mechanics, Pergamon Press (1972). 
[8] M. Fannes, J.V. Pule and A. Berbeure, Helv. Phys. Acta 55 (1982) 391. 
[9] L. Dolan and R. Jackiw, Phys. Rev. D 9, (1974) 3320. 
[10] H. Perez Rojas, Phys. Let. A (to appear) (1997). 

[11] L.D. Landau and E. M. Lifshitz, Statistical Physics, Pergamon Press (1966). 

[12] M.H. Kalos, D. Levesque and L. Verlet, Phys. Rev. A 9, (1974) 2178. 

[13] K. Huang, in Bose-Einstein Condensation, Edited by A. Griffin, D.W. Snoke and S. 
Stringari, Cambridge University Press (1995). 



8 Figure Caption 

Figure 1: (c v T)— curve obtained from our model, where the divergent behavior of the 
specific heat is shown. For T > T c the ideal gas curve has been taken. 
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